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RATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 2401 

TEMPERATURE DISTRIBUTION IN INTERNALLY HEATED WALLS OF HEAT EXCHANGERS 
WITH NONCIRCULAR FLOW PASSAGES USING COOLANTS WITH 
VERY LOW PRANDTL NUMBER 
By E. R. G. Eckert and George M. Low 


SUMMARY 

In the walls of heat exchangers that are composed of noncircular 
passages, the temperature varies in a circumferential direction. A 
prediction of the magnitude of this variation is necessary in order to 
detect the places of highest temperature and to determine the operating 
temperatures admissible for the material of which the heat exchanger is 
composed. 

A method for the determination of the temperature distribution 
within the coolant and the passage walls for a special type of heat 
exchanger is presented. The heat exchanger is composed of polygonal 
flow passages which are heated uniformly by internal heat sources. The 
Prandtl number of the coolant flowing through the passages is postulated 
to be so low that the turbulent contribution to the heat flow within the 
fluid can be neglected. The differential equation describing the heat 
flow by conduction is solved by the relaxation method. Numerical com- 
putations are carried out for a triangular and a rectangular passage 
shape. The resulting ratio of the temperature differences within the 
passage walls to the difference of the wall temperature and the fluid 
bulk temperature is considerably larger than the temperature ratio in 
similar heat exchangers using a coolant with a Prandtl number in the 
neighborhood of 1. 


INTRODUCTION 


The conventional recuperative type of heat exchanger consists of 
passages for two fluids separated by a heating surface. Heat from an 
external source is transferred continuously from one fluid to the other 
through the heating surfaces. In the regenerative type of heat 
exchanger the two fluids pass alternately through the same passages. 
During the heating period heat is transferred from a hot fluid to the 
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passage walls and is stored within the solid wall material. During the 
cooling period this stored heat is then transferred tQ a cold fluid 
passing through the heat exchanger.' 

The heat exchanger considered in this report is very similar to the 
regenerative type. It differs from it only hy the fact that the heat is 
generated internally hy heat sources located within the passage walls 
and is transferred to a coolant flowing continuously through the pas- 
sages. The heat- exchanger passages are formed hy a number of plates 
assembled to form a honeycomb. High temperatures may he anticipated 
near the corners of the flow passages thus formed. A knowledge of the 
magnitude of these temperatures is necessary in order to determine the 
admissible operating temperature of the heat exchanger. 

In a previous report (reference l) the temperature distribution in 
the walls of similar heat exchangers was studied with the assumption 
that the coolant consists of a fluid with a Prandtl number in the neigh- 
borhood of 1. In the present report, the investigation carried out at 
the HACA Lewis laboratory, is extended to include coolants with a 
Prandtl number very much less than 1. The Prandtl number is assumed 
to be so low that the turbulent contribution to the heat transport 
within the fluid can be neglected as compared with the conductive 
contribution. 

A knowledge of the velocity distribution within the passages is 
necessary in order to determine the temperature distribution. A thor- 
ough investigation of the flow of water through tubes with noncircular 
cross sections was made by Nikuradse (references 2 and 3) . The velocity 
profiles obtained in this investigation will be used in the present 
report .' 

The analysis presented applies to all passage shapes for which the 
velocity distribution is known and to all possible configurations of 
these passage shapes. Solutions are presented, however, only for heat 
exchangers composed of triangular and rectangular passages. Cross- 
sectional views of these configurations are shown in figures 1 and 2. 


SYMBOLS 

The following symbols are used in this report: 


A 


A* 


cross-sectional flow area, (sq. ft) 
A/D^, (dimensionless) 
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C 

c 

D 

d 

g 

k 

L 

L* 

IT 

n 

n* 

Pr 

P 

R 

R* 

r 

s 

T 

t B 

^w 
T * 




internal perimeter of passage, (ft) 

specific Beat, (Btu/(lb) (°F) ) 

hydraulic diameter, AA./C, (ft) 

distance between parallel plates, (ft) 

acceleration due to gravity, (ft/sec ) 

thermal conductivity of coolant, (Btu/(sec) (ft) (°F)) 

distance measured along passage wall in circumferential 
direction, (ft) 

L/D, (dimensionless) 

residual value, (dimensionless) 


coordinate normal to passage wall, (ft) 
n/D, (dimensionless) 

Prandtl number, (vpcg)/k, (dimensionless) 
summation index 
radial coordinate, (ft) 
r/d, (dimensionless) 

rate of internal heat generation, (Btu/(sec)(cu ft)) 
wall thickness, (ft) 
temperature, (°F) 


‘bulk temperature of coolant^ 





wall temperature, (°F) 
p 1 , ( dimens ionless) 

velocity components in x-, y-, and z-directions, (ft/sec) 
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u 


mean velocity. 



(ft/sec) 


u 


* 


x,y,z 


ajt Jjc £ 

x ,y*,z* 


A 

A* 


e 

v 


p 

T 

_D_ 

DT 


2 

V 


u/u, (d i mensionless) 

Cartesian coordinates 

x/d, y/D, z/D, (dimensionless) 

increment of length, (ft) 

A/D, (dimensionless) 

turbulent diffusivity of heat, (sq^ ft/sec) 
kinematic viscosity, (scl ft/sec) 
mass density, ( (lb) (sec^)/ft^) 
time, ( see) 


material derivative, 
Laplacian operator, j 



+ 


+ 



v 


<5y 


+ w 




+ 



Subscripts: 

c conditions for equivalent circular tube 

m conditions at center of flow passage 

Superscript: 

* dimensionless quantity 


ASSUMPTIONS 

A number of assumptions are necessary in order to make the problem ' 
amenable to solution and to obtain results of a general nature. The 
following assumptions are in agreement with reality for the applications 
considered: 
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1. The Prandtl number of the fluid -used as coolant is so small that 
the heat flow by the turbulent mixing movements can be neglected as com- 
pared with the heat flow by conduction. In reality the turbulent con- 
tribution to the heat flow is of some importance even in liquid metals 
used as coolants. Therefore, the calculations presented herein deal 
with a limiting case. In a previous report (reference l) analogous 
calculations were presented for a fluid with a Prandtl number of 1. 

From the results of reference 1 and the results presented herein, it 
should be possible to estimate the conditions for any fluid with a 
small Prandtl number. 

2. The passages of the heat exchanger are long enough and the Reyn- 
olds number is sufficiently high so that in the cross section investi- 
gated the flow is fully developed and turbulent. This means that the 
velocity profile does not change its shape in the direction of the tube 
axis. 


3. The property values for the coolant are constant. Under this 
condition the flow is independent of the Prandtl number and has the same 
characteristics as determined by Nikuradse for water. 

4. The rate of heat generation in the walls of the heat exchanger 
is uniform. Consequently the temperature within the coolant and the 
walls increases linearly in a downstream direction provided the flow 
is thermally developed. By extrapolation of results for fluids with 
Prandtl numbers greater than 1 (reference 4), it is found that for a 
fluid with a Prandtl number less than 1 the tube length required for 
thermal development is not greater than the tube length needed for 
velocity development. 

5. The thickness of the passage walls is small as compared with 
the hydraulic diameter. Consequently the heat conduction within the 
walls of the heat exchanger will be neglected as compared with the heat 
conduction within the fluid. Coolants with a low Prandtl number are 
fluids with a high heat conductivity value. Therefore, in al l , prac- 
tical cases the heat conductivity of the wall material is not larger 
than the heat conductivity of the fluid and the conduction in the wall 
can be neglected as long as the wall thickness is not too large. 


DIFFERENTIAL EQUATIONS DESCRIBING BEAT FLOW W ITHIN COOLANT 


The flow field within the coolant flowing through the passages is 
described by the continuity and the momentum equations. The tempera- 
ture field which is super imposed on this field is given by the energy 
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equation. As long as the property values of the coolant are independent 
of temperature and as long as free -convection effects can he neglected, 
the differential equations de scribing the flow field are independent of 
the energy equation. The flow field therefore remains unchanged when a 
heat flow from the passage walls is added to the fluid and is also inde- 
pendent of the Prandtl number. Measurements of the velocity profiles 
within noncircular passages which were made by NikuracLse (references 2 
and 3) on water therefore describe the flow field for any other fluid. 
These profiles will be used in this report. 


The temperature field within a turbulently flowing coolant is 
described by the following equation when the temperature rise by inter- 
nal friction can be neglected: 


— ta 

dt 


(— + 
\pcg 


( 1 ) 


The temperatures T in this equation are tim^e mean values and all the 
contributions of the turbulent mixing processes to the heat flow within 
the coolant are summarily contained in the term ey 2 T. The first term 


within the bracket can expressed as — by definition of the 

Prandtl number. It is known that there exists, if any, only a minor 
dependence of the turbulent diffusivity e on the Prandtl number. 
Therefore, with decreasing Prandtl number the first term in the bracket 
of equation (l) becomes more and. more important and finally for very 
small Prandtl numbers the second term eV^T can be neglected as com- 
pared with the first term. Thus the differential equation (l) reduces 
to 


2* = JL v 2 T 
DT Pr 


( 2 ) 


A Cartesian coordinate system may be arranged within the flow passage 
in such a way that the x-axis is parallel to the axis of the passage 
and normal to the plane of the drawing in figures 1 and 2. The y- and 
z-axes therefore lie within a cross section of the passage. In fully 
developed flow only velocity components in the x-direction are present. 
(The velocities connected with the secondary flow detected by Nikuradse 
and described in reference 1 are small as compared with the axial 


velocities.) Under this condition the derivative 


ST 


St 


DT 

Dt 


is, for steady 


state, equal to u The gradient ^ is a constant according to 

assumption 4. Therefore, equation (2) may be written as 
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Pr St 
~ u = 




(3) 


This differential equation describes the temperature field -within the 
fluid in any cross section of the passages. 

St 

The value of the temperature gradient ^ in equation (3) can be 

connected with the internal heat generation by a heat balance made up 
for the entire perimeter enclosing one of the fluid passages. For 
reasons of symmetry, half of the total. heat generated within the wall 
which surrounds one passage is transferred by conduction to the fluid 
flowing through this passage and increases its temperature in the axial 
direction: 

rsC —a St 

— “ PS cuA (4) 

St 

This equation is used to replace the gradient in 

S 2 T S% _ u rsC 
Sy 2 + Sz 2 u 2iA 

The boundary condition is determined by the fact that 
generated by internal heat sources within the passage 
steady state be transferred into the coolant. Since, 
assumption 5, the heat conduction within the passage walls can be neg- 
lected, this condition applies to any differential length dC along the 
circumference of i^he passages (fig. l) . Each passage wall borders on 
two fluid passages (a and b in fig. l) and the heat flow transferred 
by conduction into the coolant on both sides must be equal to the inter- 
nal heat generation within the wall. This relation is described by the 
equation 



equation (3) . 


(5) 


all the heat 
walls must in 
according to 


where n denotes the outward normal to the passage -wall surface and 
the two indices a and b refer to the two fluid passages. The inter- 
nal heat generation r per unit volume is postulated to be constant. 
Therefore, 



3?S 

-=— = constant 
k 


( 7 ) 
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Equations (5) and (7) describe the temperature field within the coolant. 


DB4ENSI0HLESS VAEUABIES 


In order to generalize the analysis,, by reducing the number of 
parameters, the equations derived in the previous section are trans- 
formed with the aid of the following dimensionless variables: 





where the hydraulic diameter D is defined as four times the cross- 
sectional flow area divided by the circumference of the passage. Equa- 
tions (5) and (7) become 


and 



( 8 ) 

(9) 


Equations (8) and (9) together with a symmetry condition fully 
determine the problem. A typical symmetry condition is shown in fig- 
ure 1. Here the normal derivative at point 1 and directed towards pas- 
sage b is equal- to the normal derivative at a similarly located 
point 2 and directed towards passage a. 


The di mensionless velocity u* is a function of y*, z*, and 

the Reynolds number and expresses the shape of the velocity field in 
the passage. For noncircular passages the velocity field is known only 
from Nikuradse's measurements for Reynolds numbers of about 40,000. 
Experiments with circular tubes have shown, however, that the shape of 
a turbulent velocity profile changes only slightly with Reynolds number. 
In addition, the temperature field is not influenced strongly by the 
velocity field (appendix A) . The dependence of u* on the Reynolds 
number will therefore be neglected, and the velocity profile will be 
regarded as a unique function for any specific passage sha|pe. Equa- 
tions (8) and (9) then have a unique solution for any given passage 
configuration. 
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CALCULATION OF TEMPERATURE DISTRIBUTIONS 

An analytic solution of equations (8) and (9) can "be found only 
when u* is a given analytic function of y* and z* and when the 
■boundaries of the passages are relatively simple shapes. Neither of 
these conditions apply to the problem under consideration. The equa- 
tions are therefore solved numerically by means of the relaxation 
method. For this purpose they are expressed in finite - difference form. 

Consider a grid, or net of points, placed into a cross section of 
the flow passages as indicated in the following sketch: 



Adjacent points are separated by a small, but finite, distance A*. The 
derivatives of T* at a point 0 can be approximated by differences of 
the temperature at surrounding points. Thus the first derivatives are 



The second derivatives are approximated by differences of the first 
derivatives and become 
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aV (T 4 *-T q *) - (T q *-T 2 *) T 2 * + T 4 * - 2T q * 

dy* 2 A* 2 A* 2 

■*V Cii*-tE 0 *) - (^o*-^5*) V + V - CTq* 

dz* 2 A* 2 ~ A^ 2 

With l^bese values, equation (8) becomes 

T x * + T 2 * + t 3* + T 4 * " 4T 0 * - 4uq*A* 2 =* 0 (10) 

The boundary condition can be expressed in finite -difference form by 
assuming, for example, that points 1, 0, and 3 lie along a passage wall 
and points 2 and 4 in two adjacent passages. The normal derivatives 
are then expressed as differences of the temperature at points 2, 0, 
and 4. Equation (9) becomes 



or 


T 2 * + T 4 * - 2T 0 * + 2A* = 0 (11) 

In some of the later calculations the wall passes through point 0 at an 
angle of 45° with respect to the grid. The boundary condition 
is then expressed as follows: 

T-j* + T g * + T 3 * + T^* - 4 Tq* + 2//2A* =0 (12) 

Equation (10) together with boundary conditions (ll) or (12) was 
solved by the relaxation method for heat exchangers composed of rec- 
tangular and triangular passages. In the case of the triangular pas- 
sages, the solution is universal, that is, it applies to all heat 
exchangers with the same passage configuration, regardless of the phys- 
ical dim ension. In the case of the rectangular passages, however, the 
solution depends on the width-to -height ratio of the passages. Because 
only one particular solution is presented for the rectangular config- 
uration, the method of solution for this case will be described in 
detail. The essential features of the relaxation method (see, for 
example, reference 5) can be outlined as follows: 

Suppose it is desired to solve a given finite-difference equation 
over a certain area of integration. The equation is of the following 
type: 
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f(T*) = 0 (13) 

The solution of this equation must also satisfy prescribed conditions 
at the boundary of the area of integration. 

First, it is necessary to select a net of points covering the 
entire area of integration. The distance A* between net points is 
arbitrary, with the accuracy of the final solution increasing as the 
distance between points is decreased. Next, values of T* are assumed 
at each net point. In general, these assumed values of the function 
will not satisfy the difference equation, and the left-hand side of 
equation (13) is equal to some residual value N instead of zero. At 
any given net point, it is then necessary to adjust T* in order to 
make N vanish at that point. This adjustment of T* also changes 
the residuals at adjacent net points. However, if this process of 
adjustment is started at the point at which the absolute value of N 
is greatest, and is then repeated for points at which the value of the 
residual is successively less, the correct values of T* for the entire 
net eventually are obtained. 

This method is now applied to a heat exchanger composed of rectan- 
gular passages as shown in figure 2. The problem can be simplified by 
taking advantage of the symmetry of the heat exchanger. Lines AB, AE, 
and HE are lines of symmetry, and there can be no flow of heat across 
these lines. Furthermore, the temperature distribution is also cen- 
trally symmetric about point C. The temperature distribution need 
therefore only be found within the area bounded by AEBCB. An enlarged 
view of this area, showing some of the net points, is presented in 
figure 3. 

The approximate finite-difference equation for all internal points 
is obtained from equation (10) . For point g it is 

Tf* + T c * + V + T k* " 4T g* - 4u g *A* 2 = N g (14) 

Point g is the point under consideration, whereas points f, c, h, and 
k are adjacent points. For points along a line of symmetry within the 
fluid, as for example point c, the following equation applies. 

o T b * + 2Tg* + T d * - 4T C * - 4 u c *A* 2 = N c (15) 

Along wall AB, which is also a line of symmetry, equation (ll) becomes 
(for point e) 


2T f * - 2T e * + 2A* = N e 


(16) 



12 


NACA IN 2401 


The following equation is found for point t and s im ilar equations apply 
for all points along wall BCD 

T p * + T u * - 2T t * + 2A* = N t (17) 

Finally, the following equation can he derived for the comers of the 
walla B and D by first writing an equation analogous to equation (17) 
for points m, x, and q, then eliminating T m * with an equation anal- 
ogous to equation (16) for points n and m, and finally combining the 
residual values into a single value N^, 

T n * + T x * - 2Tg* + 3A* = N q (18) 

Before computing the residuals at each net point, it is necessary to 
assume temperatures and to determine velocities at all points. The 
velocity distributions u* for triangular and rectangular passages were 
obtained from references 2 and 3 and are shown in dimensionless form in 
figures 4 and 5. The temperature level within any cross section of the 
passage depends on the distance of that cross section from the passage 
entrance and all temperatures found will contain an arbitrary additive 
constant. It is therefore possible to assume arbitrarily the tempera- 
ture at one point within the configuration. An estimate of all other 
temperatures can be obtained from the circular tube solution as pre- 
sented in appendix A or the flat plate solution as presented in 
appendix B. 

Once the temperature distribution has been assumed, it remains to 
determine the residuals at all net points with equations similar to 
equations (14) to (18) . The point at which N has the largest absolute 
value is then selected and T* at that point is adjusted so that the 
residual vanishes. The effect of this adjustment on the residuals at 
adjacent points is calculated with the aid of the appropriate finite- 
difference equation. The process is then carried out for the point at 
which the next largest residual appears. Eventually, if the process 
is repeated often enough, the residuals at all net points approach . 
zero, and the final adjusted values of T* satisfy the approximate 
equations at all net points. 

The accuracy of the solution can be determined by increasing the 
number of net points. The change in the temperature distribution due 
to an increase of the number of net points is a measure of the correct- 
ness of the solution. 
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BESULTS 

The method outlined in the previous section was applied to heat 
exchangers composed of a large number of triangular and rectangular 
passages. In each case the passages were staggered in order to mini- 
mize the expected hot spots. 

Temperature distribution in triangular passages. - The triangular 
configuration is represented by figure 1. Each passage in this con- 
figuration is an isosceles right triangle. Temperature distributions 
in the walls of this heat exchanger are presented in figure 6. The 
ordinate in this figure is the dimensionless wall temperature T^* 
including an additive constant. Temperature distributions found with 
nets of 28, 49, and 91 grid points are included. From a comparison of 
these distributions, it can be seen that the maximum error in the final 
solution (91 points) is less than 5 percent of the maximum temperature 
difference . 


The additive constant in the solution can be e limin a.ted by forming 
a difference of the local temperature and the bulk temperature of the 
fluid, where the bulk temperature is given by 



uT dA 


(19) 


or, in terms of the dimensionless variables 


T- 


B 



u*T* dA* 


( 20 ) 


In order to simplify the interpretation of the results of this 
paper and for purposes of comparison with the results of reference 1, 
the temperature difference (T^-T-g*) was divided by a similar temper - 


which is found in 


ature difference for a circular tube (T^-Tg ) c , 

appendix A. This temperature -difference ratio for the triangular heat 
exchanger is presented in figure 7. The use of the dimensional temper- 
atures T instead of the nondimens ional values T* as the ordinate of 
this figure is justified when the circular tube has the same rate of 
heat generation, wall thickness, and hydraulic diameter as the noncir- 
cular .passage. Contours of equal temperature (isotherms) within the 
coolant of the triangular heat exchanger are shown in figure 8. 
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Temperature distribution in rectangular passages. - Figure 2 rep- 
resents the rectangular configuration. The -width -to -height ratio of 
each passage in this configuration is 3.5. The temperature -difference 
ratio (0! v --Tg)/(T w -Tg) c for this heat exchanger is shown in figure 9. 

The temperature distribution for this case was found with a relaxation 
net of 79 points. The accuracy of this solution, however, is at least 
as good as the accuracy of the solution for the triangular heat 
exchanger, because the dimensionless distance between net points is 
smaller for the rectangular heat exchanger. Contours of equal tem- 
perature within the coolant are shown in figure 10. 

The width-to -height ratio of each passage in the rectangular heat 
exchanger is 3.5. It is to be expected, however, that the temperature 
distribution near the corners is the same for passages with any ratios 
greater than 3.5. Proof of this supposition can be found in figure 11, 
where the temperature profile within the coolant on a normal to the 
wall at point C (fig. 2) is compared with the profile between two par- 
allel infinite plates (appendix B) . It is assumed that the distance d 
between opposite walls is the same for both solutions. The comparison 
is unsatisfactory when the temperature difference is based on the bulk 
temperature of the entire passage, because that particular bulk tem- 
perature includes much of the low velocity coolant near the corners. 

If, however, the temperature difference is based on a local bulk tem- 
perature existing at the section of the passage under consideration, 
(dashed curve) then the agreement is excellent. This agreement shows 
that the influence of the corners on the flow of heat at a station 
halfway between the corners is already very small for a width-to - 
height ratio of only 3.5. Thus, the influence of one corner on the 
opposite comer is even smaller, which means that the temperature dis- 
tribution in a corner of the rectangle with the width-to -height ratio 
of 3.5 does not differ markedly from the temperature field in the 
comer of a rectangle of any greater ratio. 

Influence of velocity field on temperature field. - A certain 
error may be introduced by the assumption that the velocity profiles 
measured by Nikuradse apply under all conditions, even in noniso- 
thermal flow. In order to estimate the magnitude of this error, tem- 
perature distributions in the circular tube were calculated for both 
a constant velocity and a l/7th power velocity profile. Figure 12 
shows that the variation in the temperature distribution caused by 
this extreme change in the velocity profiles is not large. It is 
therefore believed that a small error in the velocity profiles will 
cause a negligible error in the temperature distribution in noncir- 
cular passages as well. 
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Comparison with Prandtl number 1 solution. - A comparison of the 
results presented herein and the results of reference 1 shows that the 
ratio of the temperature differences in the walls to the difference of 
the average wall temperature and the hulk temperature in heat exchangers 
with triangular and rectangular passages is very much greater for a very 
low Prandtl number than for a Prandtl number near 1. The dimensionless 
temperature difference (l^-T-g*) c for a circular tube and a coolant 

with a very low Prandtl number is 0.146 (appendix A). The corresponding 
temperature difference for a coolant with Pr = 1 is given hy equa- 
tion (31) of reference 1: 


/•m * m *\ _ 1 k*s* 

'‘-Hr -J \B > c. ~ 2 Nu 

The dimensionless temperature differences in noncircular tubes can he 
compared for Pr < < 1 and Pr = 1 hy multiplying the temperature- 
difference ratios hy the appropriate circular tube value. The dimen- 
sional temperature differences follow from the definitions of T*. It 
has to he kept in mind, however, that these definitions differ in the 
two reports, and that T* ' as used in the present report is based on 
the thermal conductivity of the coolant, whereas in reference 1, T * 
is based on the thermal conductivity of the mil material. 

A first approximation of the temperature distribution in the walls 
of heat exchangers cooled hy fluids with a Prandtl number between 0 and 
1 might he obtained hy a linear interpolation of the results of the 
present report and the results of reference 1. 


CONCLUSIONS 

A method to calculate the temperature distributiQn in a heat 
exchanger composed of noncircular flow passages and cooled hy a fluid 
with a Prandtl number much lower than 1 has been presented. 

The differential equations describing the flow of heat were written 
in terms of dimensionless variables, so that their solution becomes a 
function only of the geometry of the heat exchanger. The equations were 
then solved for heat exchangers with triangular and rectangular pas- 
sages. In the case of triangular passages, the solution is universal 
and applies to all similar heat- exchanger conf igurations composed of 
isosceles right triangles. In the case of rectangular passages, the 
solution is a function of the width-to -height ratio of each passage. 
Results for the rectangular heat exchanger were presented only for a 
width-to -height ratio of 3.5, hut it was shown that these results can 
he generalized to include passages with any ratio greater than 3.5. 
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Temper at lire differences in tlie -walls and. -within the coolant of 
these heat exchangers -were evaluated numerically. The evaluations show 
a pronounced temperature increase near the corners of the passages. 

A comparison of the ratio of the temperature differences within the 
wall to the d i fference of the average wall temperature and the coolant 
hulk temperature of heat exchangers cooled with a fluid having a very- 
low Prandtl number and of the same temperature - dif ference ratio for 
similar heat exchangers cooled with a fluid having a Prandtl number of 1 
show's that the temperature -difference ratios for the low Prandtl number 
are much greater than those for a Prandtl number of 1. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 10, 1951. 
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APPENDIX A 


TEMEERATURE DISTRIBUTION IN CIRCULAR TUBES 


The transformation of equation (8) into a cylindrical coordinate 
system with rotational symmetry results in: 


i a* a 2 T* 
51 an * 2 = 


4u* 


(Al) 


If it is assumed that half the heat generated flows into the tube 
whereas the other half flows in an outward direction, then equation (9) 
becomes 


= 1/2 = 1 <A2) 

Equation (Al) together with boundary condition (A2) can be solved 
provided that u* is a known function of R*. A simple solution can 
be obtained by making the assumption that u* is a constant and equal 
to 1. The solution of equation (Al) satisfying (A2) is 

T* = R* 2 + C ± (A3) 

where C B is an arbitrary constant. The bulk temperature as defined 
in equation (20) becomes 


T b * = 0.125 + C x (A4) 

The temperature difference in a circular tube under the assumption of 
a constant velocity, therefore, is expressed as follows 

T* - T b * = R* 2 - 0.125 ' (A5) 

A solution more representative of the problem can be obtained by 
assuming that the velocity varies as the l/7th power of (l-R) 



or, in the dimensionless system of coordinates: 

u* = § (l-2R*) l/7 


(A6) 



18 


WACA ©I 2401 


With this assumption. the eolation of equation (Al) "becomes 



In (2R* ) + (1-ZR*) 



P=1 


(1_2B*) (p+15 / 7) 

p + 15/7 


+ C 


1 


0<lf<| (A7) 

Equation (A7) has a logarithimic singularity at if=:0. A solution of 
equation (Al) at R*=0 exists., however. This solution is: 


= - 2 § - 1 L irzkr + °i (A8) 

p=»i 

The bulb temperature is obtained by multiplying equation (A7) by 
u* and integrating each term of the equation. The final solution is 

T b * = - 0.146 + C x (A9) 

A temperature-difference profile is obtained by subtracting equa- 
tion (A9) from equations (A7) and (A8) . This profile is plotted in 
figure 12, together with a similar profile as represented by equa- 
tion (A5) . 


2154 



2164 


NACA TO 2401 


19 


APPENDIX B 


TEMPERATURE DISTRIBUTION BETWEEN TWO PARALLEL INFINITE PLATES 

The following equations describe the flow of heat between parallel 
plates, extending to infinity in the z-direction and separated by a 
distance d 


d 2 T* 

, *2 
ay 


4u 1 


dr* 

ay* 


= -1 


y*=o 

D = 2d 


(Bl) 

(B2) 

(B3) 


A solution of this system of equations can be obtained if u* is 
a known function of y*. The assumption of a l/7th power profile leads 
to the following expression 


\ v y 


1/7 


or 


u* - 2 (<ty*) 1/,? 


The temperature distribution is then found to be 

T* 


to ^ 


y* + Cn 


The bulk temperature, which is defined by equation (20), becomes 

32 


_ + 

T B = 


345 


+ Ci 


(B4) 


(B5) 


(B6) 


Thus, the temperature -difference parameter for the parallel flat plates 
is 


T* - Tb* - ^ (ty *) 15 / 7 - y* 


+ 


32 

345 


0<y*< l/4 (B7) 
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Figure 8. - Temperature of coolant within passage of triangular heat exchanger. 






Temperature difference, (T*-Tg)D/d 


28 


HACA ET 2401 



Length ratio, y/d 


Figure 11. - Temperature distribution within coolant for 
rectangular configuration at point C (fig. '2) and for two 
parallel plates. 


Temperature difference 


NACA IN 2401 


29 



Radial length ratio, R* = R/D 
Figure 12. - Temperature distribution in circular tube. 
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